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Abstract 

We experimentally demonstrate a new interferometry paradigm: a self-interfering clock. We split 
a clock into two spatially separated wave packets, and observe an interference pattern with a stable 
phase showing that the splitting was coherent, i.e., the clock was in two places simultaneously. We 
then make the clock wave packets “tick” at different rates to simulate a proper time lag. The 
entanglement between the clock’s time and its path yields “which path” information, which affects 
the visibility of the clock’s self-interference. By contrast, in standard interferometry, time cannot 
yield “which path” information. As a clock we use an atom prepared in a superposition of two spin 
states. This first proof-of-principle experiment may have far-reaching implications for the study of 
time and general relativity and their impact on fundamental quantum effects such as decoherence 
and wave packet collapse. 
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Two-slit interferometry of quanta, such as photons and electrons, figured prominently 
in the Bohr-Einstein debates on the consistency of quantum theory A fundamental 

principle emerging from those debates—intimately related to the uncertainty principle—is 
that “which path” information about the quanta passing through slits blocks their inter¬ 
ference. At the climax of the debates, Einstein claimed that a clock, emitting a photon 
at a precise time while being weighed on a spring scale to measure the change in its mass- 
energy, could evade the uncertainty principle. Yet Bohr showed that the clock’s gravitational 
redshift introduced enough uncertainty in the emission time to satisfy the uncertainty prin¬ 
ciple. Inspired by the subtle role time may play in quantum mechanics, we have now sent 
a clock through a spatial interferometer. The proof-of-principle experiment described below 
presents clock interferometry as a new tool for studying the interplay of general relativity 
[3] and quantum mechanics |4j. 

Quantum mechanics cannot fully describe a self-interfering clock in a gravitational field. 
If the paths of a clock through an interferometer have different heights, then general rel¬ 
ativity predicts that the clock must “tick” slower along the lower path. However, time in 
quantum mechanics is a global parameter, which cannot differ between paths. In standard 
interferometry (e.g. 0), a difference in height between two paths affects their relative phase 
and shifts their interference pattern; but in clock interferometry, a time differential between 
paths yields “which path” information, degrading the visibility of the interference pattern 
|6j. It follows that, while standard interferometry may probe general relativity [T lSJ, clock 
interferometry probes the interplay of general relativity and quantum mechanics. For exam¬ 
ple, loss of visibility due to a proper time lag would be evidence that gravitational effects 
contribute to decoherence and the emergence of a classical world—a world of events, such 
as measurement results—as predicted by R. Penrose ra. L. Diosi im and others. 

In our experiment, atomic clocks—atoms in superpositions of internal states—pass 
through an atomic matter-wave interferometer. We demonstrate that the visibility of inter¬ 
ference patterns produced by thousands of self-interfering clocks (atoms in a Bose-Einstein 
condensate) depends on the (simulated) proper time differential between the recombined 
wave packets of each clock. We simulate the time differential or lag by artificially making 
one clock wave packet “tick” faster than the other. While our clock is not accurate enough 
to be sensitive to special- or general-relativistic effects, it is able to demonstrate that the 
proposal of Zych et al. [SJ is sound, namely, that a differential time reading affects the visi- 
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bility of a clock self-interference pattern; specifically, the visibility equals the scalar product 
of the interfering clock states. 

In principle, any system evolving with a well defined period can be a clock. In our 
experiment, we utilize a quantum two-level system. Specifically, each clock is a 87 Rb atom 
in a superposition of two Zeeman sublevels, the rrip = 1 and mp — 2 sublevels of the F = 2 
hyperfine state. 

The general scheme of the clock interferometer is shown in Fig. 1 (for additional informa¬ 
tion see [12]). To prepare the clock in a spatial superposition of two different locations, we 
make use of the previously demonstrated Stern-Gerlach type of matter-wave interferometer 
on an atom chip, creating a coherent spatial superposition of a 87 Rb BEC [13] (about 10 4 
atoms 90 /im below the chip surface). Initially, after the application of a field gradient beam 
splitter (FGBS) and a stopping pulse which zeroes the relative velocity of the two atomic 
wave packets, the wave packets are in the same internal atomic state (|F, rn F ) = |2,2) = |2)) 
as well as in the same external momentum state. The system’s external wave function is 
thus ip(x — Xi) + xp(x — x 2 ), where x* (i=l,2) are the mean values of the position of the two 
wave packets, which have the same center-of-mass momentum. A radio-frequency (RF) 7 r /2 
pulse (Rabi frequency Qr and duration T R ) tuned to the transition from |2) to |1) = |2,1) 
forms the clock by transferring the atoms from the |2) state to the internal superposition 
state (|1) + |2))/y / 2. The pulse is applied under a strong homogeneous magnetic field (AEi 2 
~ fix25 MHz) to push the transition to |2,0) out of resonance by ~180 kHz via the non-linear 
Zeeman effect, thus forming a pure two-level system for the |1) and |2) states. 

In order to examine the coherence of the clock spatial superposition, we let the two clock 
wave packets freely expand and overlap to create spatial interference fringes, as shown in 
Fig. 2(A). As two BEC wave packets are always expected to yield fringes when they overlap, 
many experimental cycles are required in order to prove phase stability or in other words 
coherent splitting of the clock. Fig. 2(B) presents the averaged picture of 100 single shots 
taken continuously over a period of about two hours. Relative to the mean of the single¬ 
shot visibility, the contrast falls by a mere 4%, demonstrating a stable phase. The phase 
distribution in the data ra reveals that the chance that the clock splitting is not coherent 
is negligible. We have thus proven with a high level of confidence that the clock has indeed 
been in two places at the same time. 

We now show that clock time is indeed a “which path” witness. For a single-internal- 
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FIG. 1: (Color online) Experimental sequence of the clock interferometer. (A) Detailed sequence 
(not to scale): Following a coherent spatial splitting by the FGBS and a stopping pulse, the system 
consists of two wave packets in the |2) state (separated in the direction of gravity, z) with zero 
relative velocity m ■ The clock is then initialized with an RF pulse of length Tr after which the 
relative “tick” rate of the two clock wave packets may be changed by applying a magnetic field 
gradient dB/dz of duration Tq. Finally, before an image is taken (in the xz plane), the wave 
packets are allowed to expand and overlap for 8 ms. See m for a detailed description of the 
sequence. (B) Each clock wave packet shows as a one-handed clock, where the hand corresponds 
to a vector in the equatorial plane of the Bloch sphere. When the clock reading (i.e. the position 
of the clock hand) in the two clock wave packets is the same, fringe visibility is high. (C) When 
the clock reading is opposite (orthogonal), it becomes a “which path” witness, and there is no 
interference pattern. 

state interferometer, a phase difference will not change the visibility of the fringes. By 
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FIG. 2: (Color online) Clock interference: (A) A single experimental shot of a clock interfering 
with itself (z axis values are relative to the chip surface). As Tq = 0 the clock rate is approximately 
the same in the two wave packets and interference is visible. As can be seen from Fig. 3, a constant 
differential rotation of the clocks, </>o, exists even for Tq — 0 (due to a residual magnetic gradient 
in our chamber |12j). This somewhat reduces the visibility. (B) To prove the coherence of the 
clock spatial splitting, an average of 100 consecutive shots such as that in (A) is presented, with 
only a minor change of visibility (44 ±1% compared to 46 db 4% for the mean of the single-shot 
visibility [12]). (C) To prove that clock time acts as a “which path” witness, we present a single 
shot in which the differential rotation angle </>o + A ujTq equals i r. Unlike standard interferometers 
in which a phase difference does not suppress visibility, and contrary to standard split-BEC inter¬ 
ference experiments in which a single shot always exhibits significant visibility, here the visibility is 
completely suppressed. (D) Similar to (B), but where one clock wave packet has been rotated by 271- 
relative to the other so that their readings are again indistinguishable (visibility is 47 d= 1%, down 
from a single shot average of 51 ±2%). The fits are a simple combination of a sine with a Gaussian 
envelope m- Throughout this work, all data samples are from consecutive measurements without 
any post-selection or post-correction. 


contrast, the relative rotation between the two clock wave packets is expected to influence 
the interferometric visibility. In the extreme case, when the two clock states are orthogonal, 
e.g. one in the state of (|l) + |2))/v / 2 and the other in the state of (|1) —12))/v/2, the visibility 
of the clock self-interference should drop to zero. We therefore apply a magnetic gradient 
pulse (inducing a “tick” rate difference A uj [12]) of duration Tq to induce a relative angle 
of rotation between the two clock wave packets (Fig. 1). When the relative rotation angle 
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is 7 r, we observe in Fig. 2(C), in a single shot, that the visibility of the interference pattern 
drops to zero. Fig. 2(D) exhibits a revival of the single-shot visibility when the differential 
rotation angle is taken to be 2tt (where we again present an average of 100 shots to confirm 
coherence). It should be noted that the differential forces induced by the gradient are not 
strong enough to break the clock apart m 

To obtain a more general view of the effect, we present in Fig. 3 the dependence of the 
interferometer visibility on the differential rotation angle between the two clock wave packets 
over the range 0 to 47t, by varying Tq to alternate between clock indistinguishability and 
orthogonality—providing “which path” information. The blue data present the clock inter¬ 
ference pattern visibility, clearly showing oscillations (consistent with the expected period 
B2|) Comparing the latter oscillations to the visibility of a single-internal-state “no clock” 
interference (£l R T R = 0; red data) confirms that the oscillations are due to the existence 
of a clock. The single-internal-state interference data also confirm that the overall drop in 
visibility is not due to the formation of the clock. This upper bound is due to the magnetic 
gradient pulse causing imperfect overlap between the two wave packets [T2J. A lower bound 
on the visibility is due to the spatial separation of the |1) and |2) wave packets (i.e. gradual 
breakup of the clock), again due to the magnetic gradient [T2] . which results in an increase 
of the visibility as expected from two independent single-state interferometers M 

The essence of the clock is that it consists of a superposition of two levels, i.e. |1) and 
|2). In Fig. 3, we chose to work with an equal population of the |1) and |2) states upon clock 
initialization to create a proper clock, thus maximizing the visibility’s dependence on the 
differential rotation. To further prove that it is the clock reading that is responsible for the 
observed oscillations of visibility, in Fig. 4 we modulate the very formation of the clock by 
varying the clock-initiating RF pulse (T R ), so that the system preparation alternates between 
a proper clock and no clock at all. Specifically, varying T R changes the population ratio of the 
two components of the clock. When the differential rotation of the two clock wave packets 
(A ujTg) is set to 7r (i.e. orthogonal clocks), as shown by the blue data, the interferometer 
visibility oscillates as a function of the ratio of the clock states’ initial population. This is 
so because when fl R T R equals multiples of ir only one of the clock states is populated and 
the system is actually not a clock. In this case we have a standard interferometer; “clock 
orthogonality” and clock time as a “which path” witness do not exist irrespective of the fact 
that AluTg = 7r, and consequently full visibility is obtained. When a clock is formed (i.e. 


6 


0.7 



FIG. 3: (Color online) Varying the orthogonality of the two clock wave packets. To study the 
properties of clock time as a “which path” witness, we measure visibility while continuously varying 
the relative rotation of the two clock wave packets (blue). Each data point is an average of the 
single-shot visibility obtained in several experimental cycles, and the error bars are the variance in 
this sub-sample. A fit returns an oscillation constant of Acj = 0.166 d= 0.003 rad//isec, consistent 
with an independent estimate (for further details see [12]). As inferred from the single-internal-state 
“no clock” interferometer (red line) the oscillations are due to the existence of a clock. Regarding 
the upper and lower bounds, see text. We note that the maximal visibility is slightly different from 
that of Fig. 2(D) as the data here are from a different run. 

when the initial populations are similar), clock time is an effective witness, and the visibility 
drops. By contrast, when AujTq = 2t r (red data), the interferometer visibility is always high 
because the two wave packets are not orthogonal whether they are clocks or states with a 
definite mp. 

Finally, we note that we consider recent works on the so-called Compton clock interfer¬ 
ometer m and the debates that ensued (see mm and references therein) to be beyond 
the scope of this paper. 

Future work will focus on clocks with increased accuracy and stability, first in the micro- 
wave and then in the optical regime, with the aim of reaching an accuracy allowing the 
experiment to detect relativistic effects. In addition, as time is considered by some a pa- 
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FIG. 4: (Color online) Varying the preparation of the clock. To further prove that it is the clock 
reading that is responsible for the observed oscillations in visibility, here we modulate the very 
formation of the clock by varying Tr, so that the system preparation alternates between a proper 
clock and no clock at all. (A) When the imprinted relative rotation between the two clock wave 
packets is i r, whether a proper clock is formed or not has a dramatic effect (blue). By contrast, 
when the relative rotation is 2 tt, whether a proper clock is formed or not has no effect (red). The 
error bars are the standard deviation of several data points. We note that in this measurement, 
there is no significant overall drop in visibility in the range of 5 oscillations as the RF pulse only 
changes the internal population. (B) The oscillation period appearing in (A) is as expected from 
an independent measurement of the Rabi oscillations induced by Tr when the rest of the sequence 
has been eliminated 


rameter which is still far from being fully understood ra, such an interferometer may shed 
new light on a variety of related fundamental questions. 
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A self-interfering clock as a “which path” witness 

Yair Margalit, Zhifan Zhou, Shimon Machluf, Daniel Rohrlich, Yonathan Japha, Ron Folman 
Department of Physics, Ben-Gurion University of the Negev, Beer-Sheva 8f 105, Israel 

This file includes supplementary Sects. S1-S7, Figs. S1-S4, and Refs. |2T1 - 127] (including 
Refs. [23H2ZI, which propose cold-atom experiments on quantum gravity—or the effect of 
gravity—that are not directly related to our work, but included for completeness). 

SI. Methods 

Our experimental procedure is as follows. We begin by preparing a BEC of about 10 4 
87 Rb atoms in the state |F, mp) = |2, 2) in a magnetic trap located 90 /im below the chip 
surface. The trap is created by a copper structure located behind the chip. (See Fig. SI.) 
The BEC atoms are released from the trap, and 0.9 ms later a field gradient beam splitter 
(FGBS) [13] is applied to create a coherent spatial superposition of the BEC. The FGBS 
consists of one radio-frequency (RF) 7t/2 pulse (of 10 /is duration), a magnetic gradient pulse 
(4 /is), and another RF pulse (10 /is). These pulses create a superposition of |2) = 12, 2) wave 
packets having different momenta. After the FGBS, we apply a second magnetic gradient 
of 90 ^s duration to zero the relative velocity of the wave packets. Clocks are initialized 1.5 
ms after trap release by a third RF pulse of duration Tr, which creates a superposition of 
the \2) and |1) = |2,1) states in each of the wave packets; then a phase-imprint magnetic 
gradient pulse of duration Tq is applied in order to change the relative “tick” rate of the 
two clock wave packets. 

The entire sequence described above is done under a strong homogeneous magnetic bias 
field of 36.7 G in the y direction, which creates an effective two-level system via the non-linear 
Zeeman effect. This field is adiabatically turned off 3.5 ms after the clock initialization (5 ms 
after the trap release), leaving earth’s magnetic field to preserve the two-level system contin¬ 
uously. After an additional 3 ms time of flight (TOF)—thus 8 ms after the trap release—we 
image the atoms by absorption imaging and generate the pictures shown in Fig. 2. 

All three magnetic gradient pulses are generated by three parallel gold wires located on 
the chip surface (Fig. SI), which are 10 mm long, 40 /mi wide and 2 /mi thick. The wires’ 


12 




10 mm 


FIG. SI: (A) A picture of the atom chip on its mount, with the copper structure visible behind 
it. Note that its orientation in the experimental setup is face down. (B) Magnetic field strength 
below the atom chip, generated by the quadrupole field via the chip wires and the bias field B y 
via external coils. The purple dot shows the location of the trapped BEC, which has in the xz 
plane a half-width of about 9 /im. (C) Schematic diagram of the relevant chip wires. Wires are 
10 mm long, 40 fi m wide and 2 /im thick. The separation of the wires’ centers is 100 /im, and 
the direction of the current / alternates from one wire to the next. The wires, being much smaller 
than the chip, are hardly visible in (A). 

centers are separated by 100 gm, and the same current runs through them in alternating 
directions, creating a 2D quadrupole field at z = 100 gm below the atom chip. The FGBS 
phase noise is largely proportional to the magnitude of the magnetic field created during 
the gradient pulse ra- As the main source of magnetic instability is in the gradient pulse 
originating from the chip, positioning the atoms near the middle (zero) of the quadrupole 
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FIG. S2: Histograms of 100 shots used to generate the averaged pictures in Fig. 2. (A) Distribution 
of phases in Fig. 2(B), with a relative clock rotation of </>o (Tq = 0). (B) Distribution of phases 
in Fig. 2(D), with a relative clock rotation of <j )o + Ac oTq = 27r. The standard deviations of the 
distributions are a = 0.454 rad and a = 0.314 rad for (A) and (B), respectively. Both distributions 
are significantly smaller than a random distribution of phases. Also shown are respective polar 
plots of phase vs. visibility (shown as angle vs. radius). 

field created solely by the three chip wires 100 /im below the chip surface reduces the phase 
noise during the FGBS operation. The chip wire current was driven using a simple 12.5 V 
battery, and was modulated using a home-made current shutter, with ON/OFF times as 
short as 1 /rs. The total resistance of the three chip wires is 13.6 Q. yielding a current of 
11.3/13.6 A ~ 0.83 A (a small voltage drop exists in the circuit itself). 

The RF signal is generated by an Agilent 33250A waveform generator and subsequently 
amplified by a Minicircuit ZHL-3A amplifier. We generate RF pulses using a Minicircuit 
ZYSWA-2-50DR RF switch. RF radiation is transmitted through two of the copper wires 
located behind the chip (with their leads showing in Fig. SI). 

S2. Coherence and visibility: further analysis 

Fig. 2(B) shows a visibility of [44 ± 1]% (error from fit). The mean visibility of 100 
single shots is [46 ± 4] % (error from standard deviation). The reduction of a mere 4% 
[=(46-44)/46] from Fig. 2(A) to Fig. 2(B) demonstrates a notable phase stability of the 
spatial interferometer. 

Figs. S2(A-B) present the phase distributions of Figs. 2(B) and 2(D), each of 100 
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FIG. S3: An average of 65 consecutive single shots, in which the phase difference between the clocks 
is 0o + A cjTq = tv as in Fig. 2C, yielding almost zero visibility. The minimal visibility 0.06) 
that persists is consistent with the result of Fig. 3, in which a n relative rotation yields a visibility 
of 0.09. Possible causes are the limited resolution (2 /is) in setting the time Tq, temporal jitter in 
the electronics that produces Tq , instabilities in the chip current, and a slight clock breakup due 
to the magnetic gradient. 

events. The probability that these phase distributions arise from a random phase process is 
negligible, which confirms the coherence of the clock interferometry. 

Fig. S3 shows an average of 65 consecutive single shots of clock interference when the 
7r relative rotation is applied, yielding very low visibility 0.06). This minimal but non¬ 
zero visibility is due partly to our limited time resolution in setting T0, timing jitter of the 
electronics, and instabilities in the chip current. See also the discussion of clock breakup in 
Sect. S6. 

S3. Fit and interpretation of Fig. 3 

The red data in Fig. 3 show a decay of the single-internal-state interference visibility 
as a function of To- We interpret this decay as the signature of non-perfect overlap of 
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the two wave packets in space. Non-perfect wave-packet overlap can occur if the atoms 
are not located directly below the center of the three-wire structure, because of initial 
fluctuations in trapping position. When the magnetic gradients are applied, they cause 
the atoms to acquire some differential momentum along other axes in space beside the z 
(gravity) direction, inducing some mismatch in the final x and y positions and thus reducing 
the interference visibility. Clearly this effect increases with Tq . In addition, the contrast of 
the clock visibility oscillations seen in Fig. 3 (blue data) is expected to decay due to the 
spatial separation of the |2) and |1) states within each of the two clock wave packets. This 
gradual breakup of the clock is due to the different force being applied during the magnetic 
gradient pulse on each of the states within the clock. It eventually results in an increase 
of the minimal visibility as orthogonality does not affect two non-overlapping single-state 
interferometers. 

In order to understand the interference of the clock on the background of the overall decay 
of visibility with Tq, we first fit the single-internal-state interference data (red points in Fig. 
3) to the expression a • e _T<3 / T1 , where a is the amplitude and T\ is the corresponding decay 
rate. We then use these values of a and T\ in fitting the clock interference data (blue points 
in Fig. 3) to the expression V (Tq) = ae~ To ^ Tl \Jl — sin 2 [((/>o + AuT G )/2]/ cosh 2 (a 0 + T g /t 2 ) 
(see Sect. S6), where «o, the initial decay of the contrast (of the clock’s visibility oscillations), 
and t 2 , the decay rate of the same contrast, are both related to the above-mentioned gradual 
breakup of the clock [the first related to a residual gradient in our system (see Fig. S4) and 
the second to the induced magnetic gradient pulse, of duration Tq] ; A uj is the relative clock 
“tick” rate; and o is the relative rotation between the two clock wave packets, again caused 
by the residual magnetic gradient in our system. 

An independent estimate of Ao; comes from calculating the expected magnetic field dif¬ 
ferences between the two wave packet positions during the magnetic gradient. From the 
Biot-Savart law, we obtain B(z + Sz ) — B(z ) = 0.0387 G. Here Ao; is 0.17 ± 0.01 rad//rs, 
where the uncertainty in Aa; is based on a position uncertainty of about 2 /.an. This estimate 
is consistent with the fit result from Fig. 3 of Aa; = 0.166 ± 0.003 r&d/ys. 

S4. Different clock preparations in Fig. 4 

While in Fig. 3 the time interval between the clock-initiating RF pulse (T R ) and the 
magnetic gradient (Tq) is minimal (8 y s), in Fig. 4 the gradient was delayed by 100 /is in 
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order to allow Tr to reach 100 //s and enable us to investigate different clock preparations. 
Longer times make the atoms fall further away from the chip and the effect of the magnetic 
gradient produced by the chip diminishes. Consequently Tq was increased to produce the 
required phase shifts. There is no significant decay observed in Fig. 4(A) as the RF pulse 
only changes the population ratio between |2) and 11), with no detrimental effects as a 
function of Tr, as can be seen in Fig. 4(B). 

In order to measure Fig. 4(B) in the same position of the atoms as in Fig. 4(A), we 
cancel the RF pulses which are part of the FGBS and set Tq = 0, while maintaining the 
magnetic-gradient pulses (FGBS and stopping pulse). We then apply the clock-initiating RF 
pulse of duration T R , and perform a Stern-Gerlach measurement of the relative population 
between |2) and |1). The Stern-Gerlach effect is created using a current pulse on the copper 
wires located behind the chip (Fig. SI). 

55. Effect of bias field inhomogeneity on the visibility 

As noted, our bias field is in the y direction. Inhomogeneities of this field in the z 
direction, along which the two clock wave packets are separated, are expected to influence 
the relative “tick” rate as well as the clock breakup. In order to calculate their effect on 
the visibility of clock self-interference, we varied the duration of the bias field during the 
clock’s time of flight (TOF) and measured the resulting visibility. We extended the TOF 
by 10 ms (total TOF = 18 ms) in order to achieve longer clock times in the bias field. Fig. 
S4 shows that the visibility of clock self-interference depends on how long the clock is in the 
bias field. The effect is analogous to the oscillations seen in Fig. 3 (as a function of Tq), but 
here Acj is much smaller since the inhomogeneity of the bias field is small. Hence, the small 
inhomogeneities in the bias field may be parameterized by «o (Sect. S3) as causing an initial 
decay of the contrast in the visibility oscillations of the clock’s interference pattern (due to 
clock breakup). These inhomogeneities also give rise to A, as they shift the differential clock 
“tick” rate between the two clock wave packets. 

56. Model for the interference of the two clock wave packets 

To calculate the interference between two separated clock wave packets, we can model 
the wave packets as Gaussian wave functions. We begin with 

e -a( Z -d/ 2) 2 | 2 ) + e -a(z+d/ 2) 2 ( S1 ) 
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FIG. S4: Effect of inhomogeneities in the magnetic bias field on the visibility of clock self¬ 
interference. We vary the duration of the clock’s stay in the bias field. The TOF is extended 
to 18 ms in order to achieve longer times in the bias field. The decrease and subsequent increase in 
visibility confirm that small inhomogeneities in the bias field induce the initial (Tq = 0) reduction 
in visibility seen in Fig. 3, and are the cause of the initial relative phase </>o of the clock wave 
packets. As the clock wave packets are separated vertically (in the z direction), the effect is due 
to a residual gradient in the z direction. 

as the state just after FGBS ra and the stopping procedure; the state (corresponding to 
the first stage in Fig. 1) is a superposition of two Gaussian wave packets with zero relative 
momentum separated by a distance d; here we neglect overall normalization. Before the clock 
initialization, both wave packets are in the internal state |2). The parameter a is a complex 
number, and its real part 3?[a] = l/2cr 2 contains the width a of the Gaussian function. Its 
imaginary part is the coefficient of the quadratic phase due to the evolution of the wave 
packets during the preparation stage, especially the evolution due to harmonic components 
of the magnetic potentials during the FGBS gradient pulse and the stopping pulse. During 
free evolution, if t 0 > 0 this quadratic phase creates focusing and then divergence of the 
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wave packet; or if to < 0, it creates divergence of the wave packet such that its consequent 
evolution is similar to the evolution of a minimal wave packet from t 0 < 0. During free 
evolution of the wave packets, the parameter a evolves as 

a(t) = ^[ctq + ih(t - to)/m\~ l , (S2) 

where <To is the wave packet’s minimal width at time t = to . The actual width a at time 
t = 0 is related to these parameters by a = cr 0 ^/l + ( hto/mo q) 2 . 

Thus, at the end of the “relative clock rotation” stage, which involves a magnetic field 
gradient dB/dz, the state is 

e ~a(z-d/ 2) 2 _|_ e ~a{z+d/ 2) 2 

where k = (nB/2h)dB/dz, Hb is the Bohr magneton, and e lkz \l) + e 2lkz \2) is the clock state. 
Note that the gradient pulse has two effects. First, it imprints different phases e lkd / 2 and 
e -ikd /2 on ^ 1C |- wo docks centered at z = d/2 and —d/2, respectively. This relative phase is 
responsible for the “which path” information that is carried by the clock. Second, it applies 
a differential momentum kick hk and 2 hk to the two states |1) and \2) that make up each 
clock state, leading to a partial breakup of the clock after propagation. In what follows we 
examine quantitatively the effect of these two aspects of the gradient pulse on the formation 
of the interference pattern. 

Equation (|S3j) can be expressed as a sum of four terms that are equivalent up to exchanges 
d -B —d, |1) B- 12) and/or k B- 2k. The time evolution of each of the terms obeys the 
Schrodinger equation and the solution is well known: 

g— a(z±d/2) 2 +ikz _ y ^—ihk 2 t/2m^—(z±d/2—hkt/rn)‘ 2 /[2(<jQ+iTi(t—to)/rn)]+ik(z—Tikt/rn) (S4) 



To exhibit the character of the interference pattern, it is sufficient that we consider the case 
h(t — to)/m Oq (the far-field limit), in which the interference pattern is fully developed. 
By detecting the atoms with an imaging light that does not distinguish between the two 
eigenstates |1) and |2), one observes a sum of two fringe patterns with detection probability 
P(z) = | (11-0) | 2 +1 (2\ip) | 2 , where | tp) denotes the overall state. In a frame of reference moving 
with the center-of-mass momentum 3hk/2, it has the form 

^ (z + A^/2) 2 
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where the two patterns have a periodicity \ t = 2nh(t — to)/md and an envelope width 
a t = cro-^/l + ( h(t — to)/ ttzcto ) 2 . The two fringe patterns, one of state |1) and the other of 
state |2), are shifted relative to one another by a phase (f) t = kdt/(t — to) and their envelopes 
are shifted with respect to each other by A z t = hkt/m = \ t 0t/2n. The shifts in phase and in 
space are one and the same and are both caused by the differential momentum kick applied 
at the clock rotation stage. If the spatial shift A z t is much smaller than the Gaussian width 
a t , then the visibility of the joint fringe pattern depends on the relative phase of the two 


fringe patterns in Eq. (S5). If the relative phase is zero or an even multiple of 7r, then the 
two interference patterns overlap and the visibility is 1. If (J) t is an odd multiple of 7r then 
the two fringe patterns are completely out of phase and the joint visibility drops to zero. In 
the case where A z t becomes considerable relative to the Gaussian width a t . each clock wave 
packet is gradually broken into its constituents, eventually forming separate wave packets 
for the states |1) and |2). In this case the cancelation of the two fringe patterns when 0 t = tt 
is not effective any more and the visibility does not drop to zero. 


In order to examine the visibility more closely, we write Eq. (S5) in a different form 
(omitting a constant pre-factor) 


P(z) oc e z2 / a t |cosh(^A^/(jj) [1 + cos(</> t /2) cos(27T2:/A t )] 

+ sinh(,sA,Zt/of) sin(0 t /2) sin(27T2;/At)| 

This form can be further simplified to the form 

P(z) oc e~ z cosh(zAz t /cr^) [1 + V(t, z) cos[27r^/A t — < p(t, z )\, 


(S6) 


(S7) 


where 

V(t,z) = \fl - sin 2 (0 t /2)/cosh 2 (^A 2 ; t /cr t 2 ) (S8) 


and 


ip = atan[tanh(^A^/cr^) tan(^/2)]. 


(S9) 


If the breakup effect is small, namely A z t <C cr t , then cosh(oAo t /cr 1 and the visibility 
V becomes V ~ |cos(0 t /2)|, which is the scalar product of the two clock states (|1) + 
e ±i<Pt/ 2 1 2 ) )/ v^2- In this case the phase cp is also negligible. The clock breakup effect makes 


20 










the visibility increase, such that even when the clock states are completely orthogonal (i.e. 
distinguishable, when <p t is an odd multiple of n) the visibility drops to l 7 ~ tanh(A^/cr t ). In 
the intermediate range where the clock breakup effect is significant but does not completely 
break the eigenstates apart, the phase <p(t, z ) causes a ^-dependent change of the periodicity 
of the fringes. 

The effect of clock breakup follows from the fact that a constant gradient was applied in 
order to simulate the time lag between the two clock wave packets. In this case a momentum 
kick to each wave packet is inevitable, even if the two wave packets are well separated at 
the time when this clock rotation is applied. In principle, the clock breakup effect could be 
avoided if the time lag could be formed by fields (magnetic or light) which have different 
amplitudes at the position of the two clock wave packets but are constant across each of the 
wave packets. 

S7. The clock interpretation 

We have interpreted our data in terms of a self-interfering clock. Nevertheless it is 
instructive to discuss an alternate interpretation that does not refer to clocks at all. For 
the purpose of this discussion, we can simplify the interfering clock wave packets and write 
them as ip+{z) and 'ip-(z), where 

il>+(z) = ip(z)e inz/x [|1) +e^ /2 |2>] /2 

r (.:) - r(.:)c [|l) + r ' o/2 |2)] /2 , (S10) 

and ip(z ) is a (one-dimensional) localized, normalized wave function. Then ip+{z) and 'tp-(z) 
represent counter-propagating and fully overlapping clock wave packets with a relative ro¬ 
tation angle 0, and their sum corresponds to the interference of the two clock wave packets. 
The visibility V of the interference pattern depends on the relative clock angle <p: 

\ip + (z)+*p_(z)\ 2 = [lp+(z) + 1 p_(z)] j [lp+(z) + 1 p-(z)\ 

= \ip{z)\ 2 1 + I e 2 ” z/A (l + e^) + ^e“ 2 ™ z/A (l + e“^) , (Sll) 

so F = 1 for (p = 2nir, for cp = (2 n + l)n the clock states are orthogonal and the interference 
pattern is 

l ^)+^- W ! 2 = l ^)| 2 ■> ( S12 ) 
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with zero visibility. In general, V = | cos(0/2)| and, as we have seen, V = 0 only when the 
clock wave packets are orthogonal. This calculation makes explicit the trade-off between 
the visibility of the interference and the “which path” information arising from the clock 

mmm ]. 

However, we can rewrite |0 + (z) + 0_( z )| 2 by collecting terms in |1) and |2). We have 

ip+(z) + ip-(z) — ip(z) [|1) cos(7r2:/A) + 12) cos(irz/\ + 0/2)] , (S13) 


hence 


\i>+{z) + i>-(z)\ 2 = \i>(z)\ 2 cos 2 (7T^/A) + cos 2 (7T^/A + 0/2) , (S14) 

which represents two superimposed interference patterns, one from the |1) state and one 
from the |2) state, with a relative phase that yields perfect visibility when 0 = 0 (up to a 
multiple of 2n) and zero visibility when 0 = 7r (up to a multiple of 27 t), in agreement with 
Eq. |S11[ yet without any reference to a clock. 

While this calculation is mathematically valid, it hides the physics of our experiment, in 
two ways. First, nothing in this calculation refers to a clock. There are innumerable physical 
realizations of a clock, and innumerable clock characteristics, such as accuracy and mass. 
The clock in our experiment is the simplest possible—a superposition of two orthogonal 
states. A superposition of N > 2 orthogonal states could be a much more accurate clock, 
and an external (rather than internal) variable could also serve to measure time [23j. For 
any possible clock there would be a mathematical analysis analogous to the one above; for 
example, corresponding to a superposition of N orthogonal states would be N interference 
patterns that could add constructively to yield perfect visibility or destructively to produce 
a flat probability distribution. What all these analyses would miss is the fact that the system 
analyzed is a clock (and therefore must measure proper time). 

Second, this analysis fails to connect the loss of visibility to “which path” information. 


In Eq. (S14) there is zero visibility when 0 = 7r, but it arises as the sum of two interference 
patterns, each with full visibility; since neither interference pattern is consistent with “which 
path” information, how could we have guessed that their sum is consistent with “which path” 
information? 

To conclude, among mathematical descriptions, the description via self-interfering clocks 
provides the best physical understanding of our experiment. 
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